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1 Introduction 

Higher-order systems appear in many models in theoretical and mathematical physics, such as 
in the mathematical description of relativistic particles with spin, string theories, gravitation, 
Podolsky's generalization of electromagnetism and others. They also appear in some problems of 
fluid mechanics, electric networks and classical physics, and in numerical models arising from the 
geometric discretization of first-order dynamical systems that preserve their inherent geometric 
structures (see [20} |2~T] for a long but non-exhaustive list of references) . 

In these kinds of systems and, in particular, in higher-order mechanics, the dynamics have 
explicit dependence on accelerations or higher-order derivatives of the generalized coordinates of 
position. So, for Lagrangian systems, if the Lagrangian function depends on derivatives of order 
k, the corresponding Euler-Lagrange equations are of order 2k. Thus, the geometric descriptions 
of these systems use higher-order tangent and jet bundles as the main tool (see, for instance, 

UUM CE3 HBl [121113]). 

Furthermore, a generalization of the Lagrangian and Hamiltonian formalisms of first-order 
autonomous mechanical systems exists that compresses them into a single formalism: the 
Skinner-Rusk or Lagrangian-Hamiltonian unified formalism, proposed by R. Skinner and R. 
Rusk for first-order autonomous mechanical systems [23] . It was generalized to non-autonomous 
dynamical systems, control systems, first-order classical field theories and higher-order classi- 
cal field theories [21 112] . The generalization of the Skinner- Rusk unified formalism for 
higher-order mechanical systems has been developed in recent papers [61 I20|. [2Tj . 

The aim of this lecture is twofold: first to review this unified formalism for higher-order 
mechanical systems and second to state the variational principles for higher-order systems and 
derive the higher-order Euler-Lagrange and Hamilton equations using this unified framework. 
These geometric variational techniques are based on those introduced for for first-order field 
theories in [H \13\ [14} I15j . Our study is made for non-autonomous higher-order mechanical 
systems (the autonomous case can thereby be obtained by using trivial bundles and removing 
the time-dependence). 

In particular, we start by introducing some basic geometrical background in Section [21 and 
then reviewing the construction of the framework for the unified Lagrangian-Hamiltonian for- 
malism for non-autonomous higher-order systems (developed in [21j) in Section [3[ The main 
contributions of the paper begin in Section H[ where we establish the variational principle and 
use it to derive the higher-order equations for the Lagrangian-Hamiltonian unified formalism, 
which are written in several equivalent geometric ways. Then, in Section \5\ these equations 
are analyzed in detail, showing how they compress not only the dynamical evolution equations 
but also the equations of the Legendre-Ostrogradsky transformation connecting the Lagrangian 
and Hamiltonian formalism, which appear as compatibility and consistency conditions for the 
equations. Other relevant results are presented in Sections [6] and where first we recover the 
generalization to higher-order systems of the Hamilton Variational Principle of the Lagrangian 
formalism and the Hamilton- Jacobi Variational Principle of the Hamiltonian formalism, and then 
the higher-order Euler-Lagrange and the Hamilton equations. All these results are obtained in a 
straightforward way from this unified formalism. Finally, some conclusions and further research 
on these topics are discussed in Section 

All the manifolds are real, second countable and C°°. The maps and the structures are 
assumed to be C°°. Sum over repeated indices is understood. 
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2 Higher-order jet bundles over IR 

(See [23] for details on jet bundles and [21] for details and proofs on the unified formalism). 

Let E — ^> R be a fiber bundle with dim.E = n+1, and let rj G 17 1 (M) be the canonical volume 
form in R. If A; € N, the kth order jet bundle of the projection tt, J k TT, is the ((k + l)n + 1)- 
dimensional manifold of the fc-jets of sections 4> 6 r(7r). A point in J k ir is denoted by j k (f>, 
where eft G r(-7r) is any representative of the equivalence class. We have the following natural 
projections: if r ^ k, 



TT. 



k . 



J k TT 



J r TT 

f4> 



TT 



k . 



J k 7T 



E 



TT 



TTOTT k : J k TT 
j k <j) 



7T(0) 



TT k , where J°tt is canonically identified with E, and ix\ 



Idjfc,,-. Furthermore, 



Notice that ttq 

if </> G r(7r) is a section of tt, we also denote the canonical lifting of 4> to J k ir by j k 4> G r(7f fc ). 

Let t be the global coordinate in K such that r\ = dt, and (t, q A ), 1 ^ A ^ n, local coordinates 
in E adapted to the bundle structure. Then, natural coordinates in J k T\ are (t, q A , q A , . . . , q A ) = 

(t,qf), with q A = (fi A , qf = — — • Using these coordinates, the local expressions of the natural 



projections are 



dt % 



ir k r (t,qf) = (t,qf) , 7T k (t,q A ) = (t,q A ) , n k (t , q A ) = (t) . 



A section ip G T{Ti k ) is holonomic of type r, 1 ^ r ^ k, if r+1 = ir k _ r+1 o ?/;, where 
cj) = TT k o ip G r(-7r); that is, the section ?/> is the lifting of a section of tt to J fc_r+1 7r. In particular, 
a section ^ is holonomic of type 1 if j fc </> = ifi; that is, ^ is the canonical A:-jet lifting of a section 
4> G r(7r), where <j) = Ti k oip. A vector field X 6 X(J k ir) is a semispray of type r if every integral 
section of X is holonomic of type r. Throughout this paper, sections that are holonomic of type 
1 are simply called holonomic. 




In natural coordinates, the local expression of a holonomic section of type r is 
Thus, the local expression of a semispray of type r is 



) Qk-r+li Wk-r+2i ■ 



X 



d 







dt +Ql dq A 



+ 



+ Qk-r+1 



d 



dq A 



+ ^k-r+l 







k—r 



dq A 



+ ... + x£ 9 



k-r+l 



dq A 



3 Lagrangian-Hamiltonian unified formalism 

Let 7r : E — > IR be the configuration bundle of a kth order dynamical system, with dim E = n + 1. 
The higher-order extended jet-momentum bundle and the higher-order restricted jet-momentum 
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bundle are 

W = J 2k - l Ti x Jk -i w T*(J fc ~M ; W r = J 2fe_1 vr x jHl J fe -V , 
where J^V = T*( J k - 1 ir)/(it k - 1 )*T*R. 

(Observe that dimT*(J fc ~M = 2kn + 2 > 2kn + 1 = dim J 2k ~ l ix = dim J^V). 
These bundles are endowed with the canonical projections 

p 2 : W-^TV^tt) 



pi : W -> J 2 *- 1 ^ 
Pjk-i^: W -> J fe_1 7r 

P>-1, : -> J fc_1 7r 



0r : W -> R 

p£: J fc -V 

Pro : VV r — >■ R . 



In addition, the natural quotient map p: T*(J fe 1 7r) — > J fe 1 7r* induces a natural projection 
p w : W -> W r . 




If (t, gj^) are local coordinates in E 1 adapted to the bundle structure, the induced coordinates 
in all these bundles are 

J 2 *" 1 * : (t, q£,..., qti,Qt, ■ ■ ■ > = (*> «iS <?/)■ 



T^J^V) : (t, g A , . . . , gj^.p^, . . . ,p h A ~ 1 ) = (t, qf,p,p l A ) 

J*" V : (f , ■ ■ ■ , ^-i,pi, • • • A' 1 ) = (*> A)- 

W: {t,q£,...,q£_ 1 ,q£,...,q£ k _ 1 ,p,p A ,...,p k A - 1 ) = (t,qf,qf,p,p A ). 

W r : {t,q^,...,qti,qt---,Q2k-l,PA,---,P A ~ 1 ) = {^qt,qf,p A ). 

Remark: The last coordinates are the real dynamical variables and so W r is the real phase 
space in the unified formalism. 

Observe that dim W = 3kn + 2 and dim W r = 3kn + 1. 



Definition 1 A section ip £ r(ppj is holonomic of type r in W, 1 ^ r ^ 2k — 1, if the section 
ipi = pi o ip G r(7f 2fc_1 ) is holonomic of type r in J 2k ^ 1/ K. 

A vector field X £ XiW) is a semispray of type r in W, 1 ^ r ^ k, if every integral section ip 
of X is holonomic of type r in W. 
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Let 9 fc „i e Q l {T*{J k ' l -K)) and = -d9 fc _i G S7 2 (T*( J^vr)) be the canonical forms 

on T*(J k ~ 1 ir). The higher-order unified canonical forms are = p^Qk-i G O l (W) and SI = 
pffik-i S 4? 2 (W). Notice that kerfi = ker^*, and then (W, Q, P^i]) is a precosymplectic 
manifold. 

In natural coordinates, the above forms are given by 

9 = p\dqf + pdt , tt = dqf A dp\ - dp A dt , 
and ker is locally given by 

ker!i= (4-'4^)- 

Definition 2 T/ie higher-order coupling 1-form C G J7 1 (W) is i/ie p^-semibasic form defined as 
follows: for every w = (y,a q ) G W, where y G J 2fc_1 7r, g = 7r^_^ 1 (y), and a q G T*( J k ~ 1 7r); if 
4> G r(-7r) zs any representative of y, and u G T W W, then 

(C(w) | u) = (a q | (T w (j fc -V o p R ))(u)} . 

As C is a /OR-semibasic form, there is a coupling function C G C°°(>V) such that C = Cp^r] = 
Cdt. In natural coordinates the coupling function is 

c= P +A4+i- (i) 

The dynamical information is introduced by giving a kth-order Lagrangian density C G 
f2 x (J k ir), which is a 7r fc -semibasic form. So we can write C = L ■ (ir k )*7], where L G C°°(J k ir) is 
the Lagrangian function. Then we denote C = (Trjf" 1 o Pl )*C. As the La grangian density is a 
7f fc -semibasic form, then £ is a pK-semibasic 1-form, and we have that C = Lp^r/ = Ldt, where 
L=(7rf- 1 op 1 )*LGC 00 (W). 

In order to have a geometric structure in W r we define the so-called Hamiltonian submanifold 

W Q = {w€ W: C(w) =CH} 4W. 

Since C and £ are both pjg-semibasic forms, the submanifold W is defined by the constraint 
C — L = 0. In natural coordinates, bearing in mind the local expression (pQ) of C, the constraint 
function is given by p + p^Qu-i — L = 0. 

From [21j we have that: 

Proposition 1 The submanifold W W is 1-codimensional, py^ -transverse and diffeomor- 
phic to W r - 

As a consequence of this, if T: W r — > W denotes this diffeomorphism, we have an induced 
section h = j o T G T(/x>v), which is specified by giving the local Hamiltonian function H = 
-L + p^q^ G C°°(Wr); that is, we have 



h(t, <£, qfiA) = (t, qf, qf, -Hj A ) = (t, qf,qf, L - p\qt x J A ) 
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The section h is called a Hamiltonian section of /iyv, or a Hamiltonian ^-section. 




J 



2k-l 



Next, we define the forms 6 r = h*Q € f2 1 (W r ) and Sl r = /i*S7 € I? 2 (W r ), whose expressions 
in natural coordinates are 



Q r = p^dqf + (L - piyqf+Jdt ; a r = dqf A dp l A + d{p l A qt +1 - L) A dt . 



Remark: The precosymplectic Hamiltonian system (W r , Q r , (Pr)*?7) ( or (W , G , (p^)*^), 
with S7 Q = j'o^) represents the higher-order dynamical system in the Lagrangian-Hamiltonian 
unified formalism. 

4 Variational Principle for the unified formalism 

Next we establish the variational principle from which the dynamical equations for the uni- 
fied formalism are derived. Our starting point is the precosymplectic Hamiltonian system 



Let r(p^) be the set of sections of p^, that is, curves ip: R — > W r - Consider the functional 



where the convergence of the integral is assumed. 

Definition 3 (Generalized Variational Principle) The Lagrangian-Hamiltonian variational 
problem for the system (W r ,f2 r , is the search for the critical (or stationary) holonomic 

sections of the functional LH with respect to the variations of tp given by tp t = ° ip, where 
{a t } is a local one-parameter group of any compact- supported vector field Z £ X v ( p ®-\W r ) , that 
is, 



(w r ,n r ,( P r R yn). 



LH: r(^) 





The main result of the calculus of variations in this context is the following: 
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Theorem 1 The following assertions on a section ip € r(p^) are equivalent: 

1. tp is a solution to the Lagrangian-Hamiltonian variational problem. 

2. tp is a holonomic section solution to the equation 

V/ i(Y)n r = , for every Y G X(W r ) ■ 

3. tp is a holonomic section solution to the equation 

i(^)(fi r o^) =0 , 
where tp' : M — > TW r is the canonical lifting of tp to T W r . 

4- tp is an integral curve of a vector field contained in a class of p^-transverse semisprays of 
type 1, {X} C X(W r ), satisfying the equation 

i{X)ft r = 0. (2) 

(Proof) We prove the equivalence 1 -4=^ 2 following the patterns taken from [11!. For the 
proof of the other equivalences, see |21| (Theorem 1). 

Let Z £ X v ^\W r ) be a compact-supported vector field, and VcKan open set such that 
dV is a 0-dimensional manifold and /o£(supp(Z)) C V. Then, 



d_ 
~dt 



t=o 



t=0 JV 



t=o Jv 



ip*a* t Q r 



ip* I lim 



a t *e r - a, 



ip*L(Z)Q r 



i)*{i{Z)&Q r + di(Z)® r ) 

*p*(-i(z)n r + di(z)e r ) 



ip*i(z)n r + / d(ip* i(z)e r ) 
ip*i{z)n r + [ *p*i(z)@ r = 

V JdV 



ip*i(z)n r , 



as a consequence of Stoke's theorem and the assumptions made on the supports of the vertical 
vector fields. Thus, by the fundamental theorem of the variational calculus, we conclude 



d_ 

dt 



J 



tp* t e r = o 



^*i(z)n r = o 



for every compact-supported Z £ JE^r^WV). However, since the compact-supported vector 
fields generate locally the C 00 (W r )-module of vector fields in W r , it follows that the last equality 
holds for every pjj-vertical vector field Z in W r . 

Now, recall that for every point w G Iin0, we have a canonical splitting of the tangent space 
of W r at w in a p™-vertical subspace and a pm-horizontal subspace, that is, 



T w W r = V w (p r R )eT w (ImtP) 
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Thus, if Y e X{W r ), then 

Y w = (Y w - T w (^ o p r R )(Y w )) + T w (i; o p R )(Y w ) = Y% + Y$ , 

with Y% G V w {p r R ) and Y$ € T w (Imtp). Therefore 

iP* i(Y)n r = ^ i(Y v )Q r + V* i{Y^)VL r = tP* i{Y^)VL r , 

since ip* i(Y v )Q r = 0, by the conclusion in the above paragraph. Now, as Y$ £ T„,(lm0) 
for every w G Ixntp, then the vector field Y^ is tangent to limp, and hence there exists a 
vector field X G X(R) such that X is ^-related with Y^\ that is, tp*X = ^^| Im ^- Then 

ip* i(Y^)Q r = i(X)ip*n r . However, as dimIm-0 = dimR = 1 and $7 r is a 2-form, we obtain that 
ip* i(Y^)Qr = 0. Hence, we conclude that ip* i{Y)Sl r = for every Y G £(W r ). 

Taking into account the reasoning of the first paragraph, the converse is obvious since the 
condition tp* i(Y)Q r = 0, for every Y G X(W r ), holds, in particular, for every Z G X V ^(W r ). 



5 Analysis of the unified dynamical equations 



In order to complete the Lagrangian-Hamiltonian unified formalism, it is necessary to analyze 
the dynamical equations. We start by using the equations written for vector fields. Thus, 
equations ([2]), with the /^-transverse condition, are 

i(x)n r = o ■ i(x)(j%)*r,?o. 

It is usual to fix the /^-transverse condition by demanding that 

i(X)(p R yr, = I . (3) 

This selects a representative in the class {X}. We will do this in the sequel. 
The first important result is |21j : 



Proposition 2 The above equations are compatible only on the points of the following subman- 
ifold of W r 

Wi = {w e W r : (i(Z)dH)(w) = 0, for every Z G kerft r j 

= G W r : (i(Y)n r )(w) = 0, for every Y G X V ^(W r )] . 



In natural coordinates, a generic vector field X G X(W r ) is given by 

X = d + f^ + Ff" + Gl 9 



dt Jl dqf J Oqf A 3p\ ' 

and the /^-transverse condition implies / ^ 0, and in particular, using ([3]), we get / = 1. 
Therefore, the dynamical equation ^ first gives 

fc-i <9L 
p A ~ 7TA = > 
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which are the compatibility relations (constraints) defining locally Wi. Furthermore, for ?J 
I < k - 1, k ^ j ^ 2k- 1; 

f A _ A . r _ dL . r i _ dL i-1 . 



therefore 



a 4 9 ^4 5 dL d dL ,-, \ 3 



at ^ ym a<^ ^ j a g / ^ a<^ d P ° A ^ \d q f J d v \ ' 

Observe that, in a natural way, X is a semispray of type k. Nevertheless, the variational principle 
requires that X must be a semispray of type 1, thus 

d 2 ^r~t a d „ 4 d dL d f dL 5 



dt f^^dqf "^4-1 dq£dp A \dqf J dp\ 



Next we must require X to be tangent to Wi. Thus, it is necessary to impose that 
LpQClvVi = ^' ^ or ever y constraint function £ defining Wi: 



dq{ 



, :2 a/. , / dL 



IV 



, ~ ^-tern^r* ,011>v "- 

a , a 



(where c^t = tt + Si+iTrr)' Repeating this procedure {k — 1 steps), we get 
at agf 

Thus we obtain a sequence of submanifolds (which can also be obtained by applying any 
other constraint algorithm [HI 116]). 

W ^ Wi ^ . . . ^ W k = W c • 
As a consequence of the last equalities we conclude that 

Proposition 3 The submanifold Wc is the graph of a map TC: J 2k ~ l ir — > J k ^ 1 n* locally 
defined by 

/, ' _r / dL \ 



FC*t = t , TCqti = qti , FLYa 1 = J2(~ 1 Y S r 



Definition 4 The map J-C: <7 2fe—1 7r — > J k ~ l ir* is the (restricted) Legendre-Ostrogradsky map. 

A Lagrangian density C is regular if the map J-C is a local diffeomorphism. Otherwise, C is said 
to be a singular Lagrangian. If J-C is a global diffeomorphism, then C is said to be hyperregular. 

In natural coordinates, the regularity condition for C is equivalent to 

( d 2 L \ 

det a b pi a )(v) ± ' for ever y y G Jkqi ■ ( 4 ) 
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Observe that X is not necessarily tangent to Wc- Thus, imposing the tangency condition 
to the last generation of constraints defining Wc, these conditions give the following equations 
(on Wc): 

(-D' fcS., - * «_,)) ^ + B-D'4 (§) - o . 

And, as a consequence of (pP, we have: 

Proposition 4 If C is regular, then there is a unique semispray of type 1, X G 3£(W r ), tangent 
to Wc, which is a solution to the dynamical equations (on Wc)- 

If C is not regular, new constraints could appear and the algorithm continues until arriving 
(in the best cases) at a final constraint submanifold Wf Wc- 

If ip(t) = (t,qf(t),q A (t),p A (t)) is an integral section of X, the above equations lead to 

■ A A -A n4 i -fl 9L .,■ 9L ,_i 

C = Qi+i ; «Sb-i = K-i ° ^ ; p A = g^A ' pa = q^a-Pa > 

and after some calculations we reach (on the points of W^, or W/), 



• o 9L . f <9L j_ x 

These equations compress both the higher-order Euler-Lagrange and Hamilton equations, as can 
be seen in the following sections. 

Remark: It is interesting to point out that the variational principle for higher-order au- 
tonomous dynamical systems, and the corresponding dynamical equations, can be obtained as 
a particular case of these results when the Lagrangian function does not depend explicitly on 
the coordinate t. 



6 Lagrangian formalism: Generalized Hamilton Principle 

In this section we show how to recover the Lagrangian formalism for higher-order mechanical 
systems. In particular, we state the classical Hamilton Variational Principle of the Lagrangian 
formalism for higher-order systems and study its relation with the unified variational Principle. 

First, consider the diagram 
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As Wc is the graph of the restricted Legendre-Ostrogadski map, we have that the map 
pf = pi o jo. '■ y^c ~^ J 2k ~ l Tt is a diffeomorphism. Then we can define the Poincare-Cartan 1 
and 2 forms in J 2k ~ 1 7r as 

&c = (jc o (pf )- 1 )*e r . ; fi £ = -de £ = (j c o ( P f )- 1 )*^ r . 

These forms can also be introduced in several equivalent ways (see, for instance, [Tl 1141 [22| 123]). 

Remark: The triple (J 2fe_1 7r, f2£, (7f 2fc ~ 1 )*?7) is the higher-order non-autonomous Lagrangian 
system associated to (W r ,Sl r , (p$i)*T]). 

Now we establish the variational principle from which we can obtain the dynamical equations 
for the Lagrangian formalism. 

Given the Lagrangian system (J 2fc_1 7r, Q,£, (7f 2fc_1 )*r/), let r(7r) be the set of sections of tt, 
that is, curves <fi: K — > E. Consider the functional 

LiTfvr) — > R 

' '(j 2fc -V)*0£ ' 



where the convergence of the integral is assumed. 

Definition 5 (Generalized Hamilton Variational Principle) The Lagrangian variational 
problem (also called Hamilton variational problem,) for the higher-order Lagrangian system 
(J 2fc_1 7r, Q,£, ( ; it 2k ^ 1 )*rj) is the search for the critical (or stationary) sections of the functional 
L with respect to the variations of <f> given by fa = &t° 4>> where {at} is a local one-parameter 
group of any compact-supported Z € X V ^(E); that is, 



d_ 

dt 



t=o 



(i 2fc -V*)*e £ = o. 



Then, as in the above section, we have: 

Theorem 2 The following assertions on a section (ft £ r( 7r ) are equivalent: 

1. (p is a solution to the Lagrangian variational problem. 

2. i(jc = j is a solution to the equation 

ijj* c i{Y)n c = 0, for every Y € X{ J 2fc ~V) . 

3. tpc = j 2k ~ l 4> is a solution to the equation 

i^'cWco^c) = , 

where R — > TJ 2k ~ l ir is the canonical lifting of ^c t° TJ 2fe_1 7r. 

4- ipc = ,? ~ "V is an integral curve of a vector field contained in a class of Tt 211 ' 1 -transverse 
semisprays of type 1, {Xc} C X( J 2k ~ l n), satisfying 



i 



(x c )n c = o. 
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(Proof) The proof of the equivalence 1 -<=>■ 2 follows the same patterns as in Theorem [TJ For 
the proof of the other equivalences, see [21] (Theorem 3). ■ 

Section solutions to the Hamilton variational problem are recovered from section solutions 
to the Lagrangian-Hamiltonian variational problem in the unified formalism. In fact: 

Theorem 3 Let ip E r( / oj K ) be a holonomic section which is a solution to the Lagrangian- 
Hamiltonian variational problem given by the functional LH. Then, the section ip£ = p\ o ip E 
r(7f 2fc ~ 1 ) is holonomic, and its projection <p = 7r 2k ~ l oip£ E r(7r) is a solution to the Lagrangian 
variational problem given by the functional L; 

Conversely, from a holonomic section ip£ = j 2k ~ x <p E r(7f 2fc_1 ) which is a solution to the 
Lagrangian variational problem, we recover a solution ip = (ipc,ipc ° FL) t° the Lagrangian- 
Hamiltonian variational problem. 




(Proof) As ip E r(/)j K ) is holonomic, then ip£ = p\ o ip E r(7r 2fc 1 ) is a holonomic section, by 
definition. 

Now, p\ being a submersion, for every compact-supported vector field X E X v ( n2k 1 )( J 2fc_1 7r) 
there exist compact-supported vector fields Y E X v ^' (VV r ) such that p\^Y = X; that is, X and 
Y are /^-related. I n particular, this holds if X is the (2k — l)-jet lifting of a compact-supported 
7r-vertical vector field in E; that is, if we have X = j 2k ~ x Z, with Z E 3L V ^(E). We denote 
by {at} a local one-parameter group for the compact-supported vector fields Y E X v( - P ^ (W r ). 
Then, using this particular choice of /O^-related vector fields, we have 




since ip is a critical section for the Lagrangian-Hamiltonian variational problem. 

Conversely, if we have a holonomic section ipc = j 2fc-1 </> which is a solution to the Lagrangian 
variational problem, then we can construct ip = (ipc, i^c °J~£), which is a section ip : M. — > Wc C 
W r of the projection p^ (remember that, in the unified formalism, the dynamical equations have 
solutions only on the points of W/;, or in a subset of it). Then, the above reasoning also shows 
also that if ip£ is a solution to the Lagrangian variational problem, then ip is a solution to the 
Lagrangian-Hamiltonian variational problem. ■ 
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In natural coordinates, if ip is given by i/)(t) = (t, qf(t),q^(t),p l A (t)), then ipc = (p\ otp)(t) = 
(t,qf(t),q^(t)), and <j>(t) = (tt 2 ^ 1 o ip c )(t) = (i,gjf(£)) satisfies the kth-order Euler-Lagrange 
equations 



dL 



dq$ 



d dL 



dt dq^ 



k d k 8L 



dt k dqi 



Finally, as a consequence of all the above results, we have the corresponding relation between 
vector field solutions to the unified dynamical equations and those which are solutions to the 
Lagrangian equations: 

Proposition 5 Let X £ X(W r ) be a vector field tangent to Wc which is a solution to the 
equations 

i(X)V r = ; i(X)(pi)* V = l. (5) 

Then there exists a unique semispray of type k, Xc £ X(J 2fc_1 7r), which is a solution to the 
equations 

i(X c )n c = ; i(X c )^ 2k - 1 )*V = l- (6) 
In addition, if £ is a regular Lagrangian density, then Xc is a semispray of type 1. 

Conversely, if Xc € X(J 2fe_1 7r) is a semispray of type k (resp., of type 1), which is a solution 
to the equations (0|) ; then there exists a unique X G X(W r ) which is a solution to the equations 
(0j and it is a semispray of type k in W r (resp., of type 1). 

(Proof) See also [5T] (Theorem 2) for a detailed proof of this statement. ■ 



7 Hamiltonian formalism: Generalized Hamilton-Jacobi Princi- 
ple 

In this section we state the Hamiltonian variational problem (Hamilton-Jacobi Principle) for 
higher-order systems, recovering it from the unified formalism. (See |21j for the proofs and 
details on the higher-order Hamiltonian formalism). 

Consider the restricted Legendre-Ostrogradsky map TC: J 2k ~ l ir — > J k ~ l ir*. First, it can be 
proved that the following statements are equivalent: 

1. Qc nas maximal rank on J 2k ^ 1 TT. 

2. T L: J 2fc_1 7r — > J k ~ l ir* is a local diffeomorphism. 

/ d 2 L \ 

3. In natural coordinates of J k ir, det I — „ j I (y) ^ 0, for every y € J k TT. 

Kdq^dq^J 

As stated in Section [5j if these conditions are fulfilled, the Lagrangan density C is said to be 
regular, and when the restricted Legendre-Ostrogradsky map is a global diffeomorphism, then 
C is hyperregular. 

Now, let V = \m(TC) -4 T*(J fc -^) and V = Im(7X) A J^tt*. If r = ir r Jk ^ o 
it k ~ l : J k ~ 1 ir* — > M is the natural projection, we denote f Q = f o j: 'P — y R. A Lagrangian 
density C € Q l (J k n) is said to be almost-regular if: 



P.D. Prieto-Martmez, N. Roman-Roy: Higher-order Mechanics: Variational Principles. 



14 



1. V is a closed submanifold of J k ~ l ir* . 

2. TL is a submersion onto its image. 

3. For every y £ J 2fc_1 7r, the fibers J-£~ 1 (J-£(y)) are connected submanifolds of J 2fc_1 7r. 

The Hamiltonian section h G r(//yv) (introduced after Proposition [1]) induces a Hamiltonian 
section h € T(fj,) defined by 

pi oh = ho P2 

Then, if @k-i and ftk-i are the canonical 1 and 2 forms of the cotangent bundle T*(J fc_1 7r), we 
can construct the Hamilton-Cartan forms in J k ~ 1 7T* and V by making 

Q h = h*e k -i € ^(J^V) ; U h = h*tt € ^(J^tt*) 
e P =j*G /l G^ 1 (7') ; fip = /fi ft € fl 3 ^) . 

Observe that FC*Q h = Q c and 7X*f4 = £l £ . 

W 





Remark: (V,£lp,T*r]) is the higher-order non-autonomous Hamiltonian system associated 
with {W r ,n r ,(p r R )*rj). 

In what follows, we consider that the Lagrangian density C G Q 1 ^ 7r) is, at least, almost- 
regular. However, all the following results also hold for regular or hyperregular Lagrangian 
densities, replacing V by the corresponding open subset of J k ~ 1 ir*, or by J k ~ 1 ir*, respectively. 

First, we establish the variational principle from which we can obtain the dynamical equations 
for the Hamiltonian formalism, and then we show how to obtain the Hamiltonian dynamical 
equations. 

Given the Hamiltonian system (V,£l-p,f*r]), let T(f ) be the set of sections of f , that is, 
curves ip : 1R — >• T 7 . Consider the functional 

H:r(f D ) — > R 

JR 

where the convergence of the integral is assumed. 

Definition 6 (Generalized Hamilton-Jacobi Variational Principle) The Hamiltonian or 
Hamilton- Jacobi variational problem for the higher-order Hamiltonian system (V, Op, T*rj) is the 
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search for the critical (or stationary) sections of the functional H with respect to the variations 
of {p given by ipt = o~t ° f, where {at} is a local one-parameter group of any compact- supported 
Z G X V ^°\V); that is 

4 / <P*t®V = (7) 



Then, as in the above sections, we have: 

Theorem 4 The following assertions on a section ip G r(f ) are equivalent: 

1. ip is a solution to the Hamiltonian variational problem. 

2. ip is a solution to the equation 

if* i(Y)Q v = 0, for every Y G %(V) . 

3. ip is a solution to the equation 

i(<f/)(npo<p) = o , 

where ip' : R — > TV is the canonical lifting of tp to TV. 

4- <p> is an integral curve of a vector field contained in a class of f a -transverse vector fields, 
{X h } C X(V), satisfying 

i{x h )n v = o. 

(Proof) The proof of the equivalence 1 2 follows the same patterns as in Theorem [TJ For 
the proof of the other equivalences, see [21] (Theorem 5). ■ 

In addition, section solutions to the Hamilton equations are recovered from the solutions to 
the dynamical equations in the unified formalism. In fact: 

Theorem 5 Let ip G r(/?™) be a critical section for the Lagrangian- Hamiltonian variational 
problem given by the functional LH. Then, the section iph = T L ° p\ ° ip = FC Q o tp/r g T(f a ) 
is a critical section for the Hamiltonian variational problem given by the functional H. 

Conversely, from a section tp^ solution to the Hamiltonian variational problem, we recover a 
solution ip to the Lagrangian- Hamiltonian variational problem. 




(Proof) Observe that J-C Q o is a submersion, since it is a composition of submersions, and 
(FC o p"[)*Qp = (p^Y^C^Qh) = (Pi)*©£ = @r- Then, for every compact-supported vector 
field Z G X V( - T °\V), there exist compact-supported vector fields Y G 3i v( - p ^ (W r ) such that 
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(J 7 C o p\)*Y = Z; that is, Z is [TC Q o p^)-related with Y. We denote by {cr[} a local one- 
parameter group for the compact-supported vector fields Y £ X v ^\W r ). Then, using this 
particular choice of (FC Q o p^)-related vector fields, we have 



d_ 
~dt 



t=o 



dt 



t=o 



(o"t ° iph)*&p 



d_ 
~dt 



t=Q 



r h L(Z)Q v = [ ^{FC op\yij{Z)® T 
Jm 



dt 



r L(Y)Q r -- 



d 

dt 



t=0 



t=0 



since ip is a critical section for the Lagrangian-Hamiltonian variational problem. 

For the converse, following the same patterns as in the theory of singular non-autonomous 
first-order mechanical systems [8] , it can be proved that there are holonomic sections ipc '■ K — > 
j2fc-i 7r Q £ p ro j ec tion Tt 2k ~ 1 such that ^ = J-Xq o ip£, and they are solutions to the La- 
grangian dynamical equations. Then, the sections ip = (tpc^h) are solutions to the Lagrangian- 
Hamiltonian variational problem (see the proof of Theorem [3]) . ■ 



In natural coordinates, if iph is given by iph{t) = (t, (^(i),^^)), ^ i ^ k 
above equations give the classical higher-order Hamilton equations: 



1, then the 



dH 



p 1 a 



dH 

8qf 



Finally, as a consequence of all the above results, we have the corresponding relation between 
vector field solutions to the unified dynamical equations and those which are solutions to the 
Hamiltonian equations: 



Proposition 6 Let X £ X(W r ) be a vector field tangent to Wc and solution to the equations 

i(X)Q r = ; i(X)(pti*ri = l , (8) 

Then there exist vector fields Xh £ X(V), which are solutions to the equations 

i(X h )Q P = , i{X h )f* n = 1 . (9) 



Conversely, if Xh £ %{V) is a vector field which is a solution to the equations ([9]), then there 
exist vector fields X £ X(W r ) which are solutions to the equations ([8]). 



(Proof) : See also |21| (Theorem 4) for a detailed proof of this statement. ■ 

Remark: It is interesting to point out that, for almost-regular systems, if the unified dy- 
namical equations have consistent solutions on a final constraint submanifold VV/ W r , then 
the Lagrangian and Hamiltonian equations have consistent solutions on final constraint subman- 
ifolds Sf = p\(yVf) "—^ J 2fc_1 7r and Vj = V, respectively. Then all the above results 
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hold on the points of these submanifolds instead of W r , J 



2k-l 



it, and V, respectively. 




J 



■2k-l 




8 Conclusions and further research 

We have made an accurate revision of the generalization of the Lagrangian-Hamiltonian unified 
formalism of R. Skinner and R. Rusk to higher-order dynamical systems. We have analyzed the 
non- autonomous case, since the autonomous case can be obtained as a particular situation of this. 
This particular situation consists in using trivial bundles and removing the time-dependence (see 
[20]). This unified formalism constitutes a nice framework which allows us to study different 
kinds of problems in a simpler way. In particular, singular (constrained) systems can be analyzed 
more easily. 

In particular, as a new contribution, the classical variational principles of first-order mechan- 
ics are generalized to this framework, in order to state the dynamical equations for higher-order 
mechanics in several equivalent ways. 

Therefore, the Lagrangian and Hamiltonian structures, equations and solutions of higher- 
order mechanics are recovered from those obtained in the unified formalism, which also includes 
the corresponding Lagrangian and Hamiltonian variational principles: the generalized Hamilton 
and Hamilton-Jacobi Principles respectively. 

Several interesting physical examples have been studied using this formalism; for instance 
the Pais-Uhlenbeck oscillator and the shape of a deformed elastic cylindrical beam with fixed 
ends, as regular systems; the second-order relativistic particle, first as a free particle and later 
subjected to a potential, as singular systems [201 EI] , and also underactuated control systems 



This generalization of the Lagrangian-Hamiltonian unified formalism to higher-order dynam- 
ical systems using a general fibre bundle E over R (instead of the classical approach using trivial 
bundles) is a first step towards the study of higher-order classical field theory. However, replacing 
the base manifold M with an orientable m-dimensional manifold M gives rise to new difficulties, 
such as defining a suitable fiber bundle to act as the phase space; obtaining the Legendre map 
without ambiguities, or obtaining the relation between the momenta (which is crucial in our for- 
mulation). Nevertheless-I-, our future aim is to obtain an unambiguous Lagrangian-Hamiltonian 
unified formalism for higher-order classical field theory, thus completing previous works [H [25] . 
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